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EFFECT OF SATELLITE SPIN ON EXPLORER 33 
AND 35 DOPPLER TRACKING DATA 

J. W. Marini 
C. W. Murrary, Jr. 


ABSTRACT 

Tracking data obtained by the Goddard Range and Range Rate 
Tracking Stations from Explorer 33 and 35 Satellites is analyzed. 

The residuals of the range rate data contain periodicities related 
to the spin of the satellites and the data sampling rate. Analysis of 
the satellite antenna indicates that the periodicities are caused by 
variations in the phase pattern of the antenna, which can also produce 
a fixed bias in the range-rate measurement. 

Using a non-linear model the spin rate can be determined from 
the data. For the two examples given in the report, the rms of the 
residuals from the non-linear model is reduced by a factor of 1/2 
and 1/3 compared to the rms of the residuals from a least squares 
straight line fitted to the same data. 
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EFFECT OF SATELLITE SPIN ON EXPLORER 33 
AND 35 DOPPLER TRACKING DATA 


SUMMARY 

Tracking data obtained by the Goddard Range and Range Rate Tracking Sta- 
tions from Explorer 33 and 35 Satellites is analyzed. Both satellites are spin 
stabilized and employ a canted turnstile antenna. 

The range rate data are recorded at the sampling frequency of one meas- 
urement per second. Each measurement consists of the time (approximately a 
half second) required to count a fixed number of cycles of the difference between 
the received and transmitted frequencies (plus a fixed reference frequency). The 
data is consequently a measure of the total change in frequency in the round 
trip of the signal from the ground station through the satellite and back. The 
two-way Doppler shift produces a slowly-varying contribution to this total 
change in frequency. When this contribution is removed by fitting a smooth 
curve (a straight line) to the data, the residuals contain periodicities which can 
be attributed to beat frequencies between the harmonics of the satellite spin rate 
and multiples of the sampling frequency. 

Using a simple model for the satellite antenna configuration, it is shown that 
the appearance of the harmonics is caused by the presence of the same harmonics 
in the phase pattern of the antenna. The analysis also shows that the satellite 
spin can cause a constant contribution to the total frequency change measured 
which would appear as a fixed error in the Doppler frequency calculated. 

The spin rate can also be determined by fitting the data, in a least squares 
sense, with a linear term plus sinusoidal terms. The fundamental sinusoidal 
frequency extracted is the satellite spin rate. This spin rate determination 
agreed to within 0.2% of that examined from references 2, 3, 6, and 7. The 
accuracy to which the spin rate can be determined in this manner is a function 
of the number of data points, the assumed model, and the short term fluctuations 
on the data (Appendix C). By subtracting this non-linear model from the raw 
data, residuals are reduced by as much as 70% relative to that obtained by a 
linear fit. Thus an error model of this form when subtracted from the raw 
tracking data will result in a reduction in the uncertainty of the computed orbit. 
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INTRODUCTION 


Spin stabilization is one of the methods for controlling the attitude of a 
satellite, i.e., its orientation with respect to some reference frame. By this 
technique, one of the body axes is fixed in inertial space while the satellite 
rotates around it. 

Explorer 33 (AIMP-D) and Explorer 35 (AIMP-E) are both spin stabilized 
for attitude control. 

Explorer 33 was launched 1 July 1966. For the data epochs considered in 
this report, the satellite was in an elliptical orbit about the earth with the or- 
bital parameters varying rapidly due to lunar perturbation. In general, the 
orbital parameters remained within the following bounds (references 1 and 2). 


Perigee radius: 20,000 km to 130,000 km 

Apogee radius: 400,000 km to 530,000 km 


Explorer 35 was launched 19 July 1967 and is in an elliptic orbit about the 
moon with an aposelene radius (greatest distance from the center of the moon) 
varying between 9,300 km and 9,600 km and a periselene radius (shortest dis- 
tance from the center of the moon) varying between 2,500 km and 2,700 km (ref- 
erence 3). Both these satellites have a canted turnstile antenna located sym- 


_ j j _ i 


- trrnn 


meiricauy wim respect to me spin axis anu uom operate at vnr tor purpuoeo 
of obtaining range and range rate information. 


The raw data received at the Goddard Range and Range Rate Tracking 
Stations from these satellites and analyzed in this report is the time T in 
seconds required to count a fixed number of cycles, N = 14,328, of the round- 
trip frequency shift f d of the uplink carrier frequency* plus a known frequency 
offset, f b = 30 KHz (reference 5). Since f d is ordinarily much less than f b , 
the time T is ordinarily about one half second. The data are collected at the 
rate of one measurement per second. 


*The satellite transponder actually receives at one VHF frequency and transmits at another, mak- 
ing use of an onboard oscillator whose effect is removed on the ground in the processing of the 
received signal (reference 4). The overall effect of the signal processing, however, is to pro- 
duce on the ground the same (except for the way that spin alters the measurement) result that 
would be obtained by coherently turning around in the ground transmitted carrier at the same 
frequency. 
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While the frequency offset f b is constant, the frequency difference f d , which 
is composed almost entirely of the two-way Doppler shift between the ground 
station and the satellite, may change somewhat over the half second counting 
interval. Since the frequency counted over this interval is the sum f d + f b , and 
since the average number of cycles per second is, of course, N/T, the average 
value of f d is given by 


f„ = L - (N/T) - f, 


(1) 


Since f d < < f b , using (1), T can be approximated by a linear function of f d . 


T = (N/f b ) - (N/fg) f d . 


( 2 ) 


It is shown in Appendix A of this report that the phase pattern of the turn- 

Diuc tuiiuuaa uou uo CAprcDDUu ao a uiu;cu iuuuuv/u |nuo a xv/vjav/ iuuv/muh kjx 

the angle of rotation. As the satellite spins, the phase pattern introduces a 
bias in the two-way difference frequency f d and hence T from equation (2). The 
bias results from the linear term in the phase pattern. Harmonics of the spin 
rate are caused by the periodic terms. 


Since the data is sampled, the frequencies which appear are the beat fre- 
quencies between multiples of the sampling rate and harmonics of the spin rate 
f s . Due to both geometrical and electrical symmetry of the antenna, the even 
harmonics of the spin rate predominate. 


Although the sampling is not instantaneous, the periodicities still appear in 
the raw data. The counting of N cycles amounts to an integration of the instan- 
taneous difference frequency, and if the latter is expanded in a Fourier series 
and integrated over the sampling period, term by term, the new series generated 
by this procedure contains exactly the same harmonic terms, although with dif- 
ferent amplitudes and phases. 
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ANALYSIS 


Figuro 1 is a typical plot of tho raw data analyzed in this report, i.e., the 
time in seconds T required to count a fixed number of cycles N (or zero cross- 
ings) of tho two-way difference frequency f d plus a reference frequency f b . 
From inspection of the figuro, it can bo seen that there is a linear trend in the 
data. 


The method of analysis used in this report (described in Appendix B) con- 
sisted in fitting a least squares straight lino to the raw data, then determining 
the frequency content in the data by spectral analysis of the residuals from 
this straight line. Figure 2 is tho spectral density of the residuals from a 
least squares straight lino fitted to the raw data plotted in Figuro 1. 

It should be noted that a least squares straight line provides a better fit to 
the raw data analyzed in this report than a quadratic or higher order polynomial 
in the sense that the root mean square of tho residuals is not improved by In- 
cluding higher order terms and that spectrum analysis of the residuals from a 
higher ordered polynomial (through 6th order) reveals the same periodicities. 

It can be seen from inspection of Figure 2 that the following frequencies 
appear in the spectrum: 0.1730, 0.3080, 0.3460, 0.4135, and 0.4810 cycles/ 
second. 

From reference 6, the spin rate of Explorer 35 for 18 February 1968 was 
obtained as 24.81 rpm, or 0.4135 cycles/second. The peak in the spectrum at 
0.4150 cycles/second, therefore, is the fundamental of the spin frequency. The 
appearance of the other frequencies in the data is due to the phenomenon of 
aliasing (Appendix B). 

Since the sampling rate is l/second, the maximum resolvable frequency is 
0.5 cycles/second. Therefore, the frequencies 0.1730, 0.3080, 0.3460, and 0.4810 
cycles/second which also appear in the spectrum are the beat frequencies be- 
tween multiples of the sampling rate and harmonics of the spin rate f s . 

Twice the spin rate* beating with one multiple of the sampling rate is 


1 - 2 f g = 0.1730 eye 1 es/ second. (3) 

which corresponds with 0.1725 cycles/second. 

*Here we ore using for f the spin rate obtained from reference 6. 


4 





EXPLORER 35 ROSMAN, N.C. 
EPOCH 18 FEBRUARY 1968 






Four times the spin rate beating with two multiples of the sampling rate is 


2 - 4f fl s 0.3460 cycles/second, (4) 

which corresponds with 0.3425 cycles/second. 

Six times the spin rate beating with two multiples of the sampling rate is 


6f s - 2 = 0.4810 cycles/second, (5) 

which corresponds with 0.4875 cycles ''second. 

Finally, eight times the spin rate beating with three multiples of the sampling 
rate is 


8k s - 3 = 0.3080 cycles/second, 


(6) 


which corresponds with 0.3150 cycles/second. 

Table I contains a Summary of the Data Analysis, giving satellite identifica- 
tion, epoch of the data (in Greenwich Mean Time), Goddard Range and Range 
Rate Tracking Station, number of points :n the data stretch, spin rate and spin 
period of the satellite obtained from references 2, 3, 6, and 7, beat frequencies 
between multiples of the sampling rate and harmonics of the spin rate, and 
corresponding frequencies appearing in the spectrum of the residuals from a 
least squares straight line fitted to the data. 

From inspection of Table I it can be seen that the even harmonics of the 
spin rate predominate. As shown in Appendix A, this is due to both geometrical 
and electrical symmetry of the antenna configuration used on Explorers 33 and 
35. 


APPLICATION OF A NON-IINEAR MODEL TO THE DATA 

Having determined through spectral analysis which beat frequencies pre- 
dominate, this information can be used to construct a more accurate model for 
the data. By applying a nonlinear regression technique (reference 8), least 
squares estimates for the parameters of the model can then be obtained. There 
are two results: 
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1. The spin rate of the satellite can be determined from the data in a least 
squares sense. 


2. The rms of the residuals from this non-linear model will be reduced 
over that of a straight line. 

In order to demonstrate the above, the following non-linear model was 
applied to raw data taken at Rosman, N.C., on 18 February 1968, 6 hours, 36 
minutes, 29 seconds GMT, and at Carnarvon, Australia, on 20 February 1968, 
17 hours, 17 minutes, 34 seconds, GMT, from Explorer 35 (AIMP-E). 


T, = Cj + c 2 i + ajSin^M) + b 1 cos(27rf g i) + a 2 sin(27r(l - 2f g )i) 


+ b 2 cos(2tt(1 


- 2f g )i) + a 3 sin(27r(2 - 4f g )i) + b 3 cos(27r(2 - 4f s )i) 


+ a 4 s in(277(6 f g - 2)i) + b^cos(2rr(6 i 3 - 2)i) + a s sin(2^(8f g - 3)i) 


(7) 


+ b 5 cos(27r(8 f g - 3)i) 


(i = 1, 2, • • • , 300) 


where Cj , c 2 , a t , a 2 , a 3 , a 4 , a g , bj , b 2 , b 3 , b 4 , b g , f g are the parameters 
of the model. 

For each data stretch 300 points were available. A generalized least 
squares program (reference 8) was used to find the best estimates of the 
parameters in the above model. That is, the parameters were solved for in 
the program to minimize the sum of the squares of the residuals of the data 
points from the model. Results of the analysis can be seen in Table II which 
shows the solution for each parameter with its associated uncertainty or 
standard deviation. (See Appendix C for the standard deviation of a least 
squares estimate in a non-linear regression model.) 

From Table II the solution for the spin rate of Explorer 35 on 18 February 
1968 using the Rosman data is f 8 = 0.41440654 cycles/second (24.8644 rpm) with 
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Table II 

Comparison of Non-Linear Model With Straight Line 
Fitted to Explorer 35 Data 
(Least Squares) 


Station: Rosman, N.C. 

Epoch : 18 February 1968, 6 hours, 36 minutes, 29 seconds GMT. 
Number of Points : 300 
Straight Line Model : 


T, = a 0 + a,i 

(i = 1, 2, 3 300) 


Standard Deviation of 
Estimate of Parameter: Parameter Estimate 

a n 0.47579445 sec o* (0.57549889) 10" 6 sec 

U Bq 

a t (0.13103713) 10" 5 (dimensionless) cr~ (0.33143597) 10 -8 (dimensionless) 
RMS of residuals = 0.49715118 x 10 " 5 seconds (32 cm/second) 


Non-Linear Model: 


Tj = Cj + c 2 i + a 1 sin(277f i) + bjCos(27rf i) + a 2 sin(27r(l - 2f g )i) 
+ b 2 cos(277(l - 2f g )i) + a 3 sin(27r(2 - 4f g )i ) 

+ b 3 cos(277(2 - 4f s )i) + a 4 sin(27r(6f s - 2)i) 

+ b 4 cos(27r(6f g - 2)i) + a 5 s in(27r(8 f s - 3)i) 

+ b 5 cos(27r(8f s - 3)i) (i = 1, 2, 3, • ■ • , 300) 
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Table Il-(Continued) 


a 


Estimate of Parameter: 


Standard Deviation of 
Parameter Estimate 


c x 0.47579569 sec 

c 2 (0.13106123) 10" 5 (dimensionless) 
aj (0.68193717) lO’ 6 see 
f 0.41440654 (cycles/sec) 
b, (0.81426620) 10" 6 sec 
a 2 (0.20113182) 10" 5 sec 
b 2 (0.40956804) 10- s sec 
a 3 -(0.16844185) 10" 5 sec 
b 3 (0.25006524) 10" 5 sec 
a 4 (0.12792709) 10' 5 sec 
b 4 (0.12415350) 10' s sec 
a 5 -(0.15831928) 10' 5 sec 
b 5 (0.10895869) 10~ s sec 


cr* 

c , 


o> 

c, 


a- 

n, 


. 




CT« 

n. 


err 


C7~ 

n. 


CT 6. 


cr-' 

a. 


a* 

n. 


(0.27217461) 10'° sec 
(0.15753571) 10 " 8 (dimensionless) 
(0.19377216) 10" 6 sec 
(0. 1503465^) 10- 4 (cycles/sec) 
(0.19271141) 10" 6 sec 
(0.22455975) 10" 6 sec 
(0.20092015) 10 ~ 6 sec 
(0.24293704) 10 “ 6 sec 
(0.21283401) 10 " 6 sec 
(0.21874661) 10' 6 sec 
(0.21941388) 10" 6 sec 
(0.22593718) 10 ' 6 sec 
(0.264495,28) 10 ' 6 sec 


RMS of residuals = (0.23628070) 10' 5 seconds (15 cm/sec). 


Station : Carnarvon, Australia 

Epoch: 20 February 1968, 17 hours, 17 minutes, 34 seconds GMT. 
Number of Points: 300 
Straight Line Model: 


T. 


a. 


V 


(i = 1, 2 , * • • i 300) 
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Table Il-(Continued) 


Standard Deviation of 
Estimate of Parameter: Parameter Estimate 

. 0,48093082 sec (0.77376408) 10 " 6 sec 

u n 0 

j (0.90051096) 10" 6 (dimensionless) a* (0.44561901) 10" 8 (dimensionless) 

i 

IMS of Residuals = (0.66842480) 10' 5 seconds (42 cm/sec) 

Non-Linear Model: 


4 r Cj + c 2 i + aj sin(27rf s i) + b 1 cos(277f g i) + a 2 sin(27r(l - 2f g )i) 

+ b 2 cos(27;(l - 2f g )i) + a 3 sin(27r(2 - 4f g )i) + b 3 cos(27r(2 - 4f g )i) 
+ a 4 sin(27r(6f g - 2)i) + b 4 cos(2 7 r(6 f g - 2)i) + a 5 s i n(27r(8 f g - 3)i) 
+ b 5 cos(27r(8 f g - 3) i) 


Estimate of Parameter: 

j (0.48093177) sec 
2 (0.90050782) 10" 6 (dimensionless) 
! -(0.61840445) 10 " 6 sec 
f g (0.41435270) cycles/sec 
bj -(0.54222824) 10“ 6 sec 
a 2 -(0.60332663) 10" 5 sec 
b 2 -(0.40009164) 10 " s sec 
a 3 (0.47320679) 10" 6 sec 


(i = 1, 2, 3. . . . , 300) 

Standard Deviation of 
Parameter Estimate 

cr s (0.23114225) 10" 6 sec 

(0.13378799) 10 ” 8 (dimensionless) 

2 

o-- (0.16438325) 10' 6 sec 

(0.98549050) 10' 5 (cyeles/sec) 

8 

o-g (0.16351419) 10' 6 sec 

o* (0.17983013) 10" 6 sec 

ot (0.19779028) 10" 6 sec 
2 

cr s (0.22675795) 10 ” 6 sec 
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Table Il-(Continued) 


Estimate of Parameter: 


Standard Deviation of 
Parameter Estimate 


b 3 (0.41969521) IQ' 5 £ 
a 4 -(0.10093327) 10 " 5 
b 4 -(0.245236GG) 10-5 
a 5 -(0.13510547) 10" 5 

b g -(0.1080719G) 10- 5 

RMS of Residuals = (0 


ec 

<7p 

“3 

sec 

<7- 

n 4 

sec 

(7f 

b 4 

sec 

a- 

n 5 

sec 

% 

20065898) lO"* 5 i 

jecor 


(0.1G472243) 10” 6 sec 
(0.21249993) 10“ 6 sec 
(0.17110756) 10 " 6 sec 
(0.18353051) 10 “ 6 sec 

(0.19343074) 10" 6 sec 

is (13 cm/sec) 


an uncertainty of 0.1503 x 10' 4 cycles/second (0.0009 rpm). The solution for 
the spin rate of Explorer 35 on 20 February 1968 using the Carnarvon data is 
f 8 = 0.41435270 cycles/second (24.8612 rpm) with an uncertainty of 0.0985 x 10" 4 
cycles/second (0.0006 rpm). 

Also, the rms of the residuals from the model using the Rosman data is 
0.23628070 x 10“ 5 seconds (or 15 cm/sec) while the rms of the residuals from 
a straight line fitted to the same data is 0,49715118 x 10“ 5 seconds (or 32 cm/ 
sec). The rms of the residuals from the model using the Carnarvon data is 
0.20065898 x 10" 5 seconds (or 13 cm/sec), while the rms of the residuals from 
a straight line fitted to the same data is 0=66842480 x 10" 5 seconds (or 42 cm/ 
sec). 

Application of a slightly different model to account for only the second and 
fourth harmonics of the spin rate had little effect on the least squares solution 
for the spin rate. However, the rms of the residuals from this latter model 
was higher than the rms of the residuals from the model in (7). 


CONCLUSIONS 

Results of a spectrum analysis of the residuals from a least squares 
straight line fitted to two-way Doppler data (the time in seconds required to 
count a fixed number of cycles of the two-way difference frequency plus a ref- 
erence frequency) received at Goddard Range and Range Rate Tracking Stations 
from the Explorer 33 (AIMP-D) and Explorer 35 (AIMP-E) satellites reveal that: 

1. The spin rate and harmonics of the spin rate can be detected in the data. 
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Using a simple model for the VHF turnstile antenna configuration on these 
satellites, it is shown that tl j antenna phase pattern can be expressed as a linear 
function plus a periodic function of the angle of rotation. As the satellite spins, 
the varying phase pattern introduces a bias in the two-way difference frequency. 
Results of the antenna analysis indicate that: 

2. A bias of several centimeters per second in the two-way Doppler data 
results from the linear term in the antenna phase pattern. 

3. Harmonics of the spin rate are caused by the periodic terms in the 
antenna phase pattern. 

Applying a non-linear model to the raw data, the spin rate of the satellite 
can be determined in a least squares sense. For the two examples given in the 
report, the rms of the residuals from the non-linear model was reduced by a 
factor of 1/2 and 1/3 compared to the rms of the residuals from a least squares 
straight line fitted to the same data - from 32 cm/sec to 15 cm/sec and from 
42 cm/sec to 13 cm/sec. 
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APPENDIX A 


Periodic Properties of a Rotating Turnstile Antenna 


Both Explorer 33 and Explorer 35 satellites make use of a canted turnstile 
antenna (Reference 9) consisting of four quarter-wave monopoles spaced sym- 
metrically about the spin axis of the satellite as shown in Figure Al. The 
monopoles are fed with relative phase angles of 0°, 90°, 180° and 270° so that 
opposing terminals are 180° out of phase. The exact radiation pattern of such 
an antenna would be difficult to calculate in view of the irregular shape of the 
body of the satellite and the presence of the booms and the solar paddles. Some 
general properties of the radiation pattern, nevertheless, can be derived by 
examining a simple model of the antenna system. 

The model adapted for this purpose is a pair of crossed center-fed dipoles, 
driven in phase quadrature. Such an arrangement preserves the phase relation- 
ships used in the satellite antennas , since a center-fed dipole is equivalent to 
two opposing monopoles fed 180° out of phase. The dipoles will be assumed to 
be no longer than a half wave in length. 

The model is illustrated in Figure A-2. The dipoles, each of length 2 4, lie 
in the x-y plane. 

The currents on the dipoles, l x (x) and I 2 (y), are assumed to be standing 
waves. They can be written as 


Ij(x) = Ij [l x (x)/l t ] a Ijf(x) 
I 2 (y) = I 2 U 2 (x)/I 2 ] - I 2 f(y) 


(Al) 


where the normalized functions f (x ) and f (y ) describing the shape of the 
currents on the dipoles have identical forms since the dipoles are assumed to 
have equal lengths. l x and I 2 are defined to be the currents at the center of 
the first and second dipoles respectively so that 


f(0) = 1 


f(4) = 0. 


(A2) 
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Because of the geometrical symmetry it will also be true that f is an even 
function 


f(x) - f(-x). 


(A3) 


The customary assumption that f(x) has a sinusoidal form is not needed in the 
analysis that follows. 

The usual spherical polar coordinates r , 6 and 0 are used to describe the 
electric intensity of the radiation field of the antenna system. This electric 
intensity will have the components E 0 ( r» 6><P) and ( r, $,<P) only, since E f ( r , 6,<P) 
is negligible in the far field. 

The field components are easily derived using standard procedures 
(Reference 10). The components are 


e ~jkr 

E* =- j7?k cos 6 

2ttt 


I 2 sin 0 I 

•'O 


f (y)cos(ky s in# s in0)dy 


(A4) 


f 


+ IjCos0 I f(x) cos(kx s in# cos 0)dx 


e <£ = -j T?k 


jkr 

2rr r 


I 


+ I 2 cos 0 I f (y) cos (ky sin 8 sin0)dy 


- Ij s in 0 I 
•'o 


(A5) 


f (x) cos (kx sin #cos 0) dx 


in which use has been made of (A3). Here r) is the characteristic impedence of 
free space and k is the wave number k = 2 tt/K where k is the wavelength of the 
emitted radiation. The factor exp j cot is understood. 

Assume that the current at the terminals of the second antenna leads that at 
the terminals of the first antenna by 90° so that the antennas are fed in 
quadrature; and assume also that the magnitude of the feed currents are equal. 
Then 



(A6) 
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Equations A4 and A5 become 


>- jkr 


E« = 1. 97 k COS 6 

2 7TT 


in0 J 


sin0 l f(y)cos(ky sin 0 sin 0) dy 
0 


1 


(A7) 


- j cos 0 I f(x) cos(kx sin 0 cos 0)dx 


E 0 = Ij^k 


B - jkr 

2tty 


( 

j s in 0 J* 


cos 0 | f(y) cos (ky sin 6 sin 0)dy 
"0 

+ j sin 0 | f(x) cos(kx sin 6 cos 0)dx 
0 


(A8) 


The purpose of this analysis is to examine the behavior of the phase pattern 
of the antenna as a function of the angle 0 , since, as the antenna spins, the field 
point relative to the antenna has fixed values of r and 6 , but a varying value of 
0. By direct substitution into equations (A7) and (A8) it can be seen that 


l* 

77\ 

= j E 9 (0) =E B (0)e i ” /2 

(A9) 

\ 


{* 

7 T\ 
+ 2) 

II 

II 

S 

'S 

\ 

to 

(A10) 


Fixing attention on E^ (similar arguments apply to the component), and 
writing this component in polar form 


E 0 = A «.)ei «<*>. 


(All) 


Equation (A10) requires that 


A (0 +-— ) = A (0) 

a ^0 + ~- j = a (0) + + 27rn 


(A12) 
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where n is an integer. In the case under consideration, where the dipoles are a 
half wave length or less, n must be equal to zero. This follows because both of 
the integrals appearing in (A8) are positive for all values of 0 since none of the 
factors in the integrands goes negative over the range of integration. A non-zero 
value for n would require that the complex number in the brackets in (A8) must 
circle the origin as 0 is increased by v/2 , but this is impossible since cos 0 + 
j sin 0 does not circle the origin over this increment. 

Consequently 

a + t) = + y ( A14 ) 


whence a(0) -0 is periodic with period rr / 2 and therefore by expanding in 
Fourier series, a(0) may be written in the form 


a(0) = a. 0 + 0 + a t sin(40 + /3 X ) + a 2 sin(80 + fS 2 ) + . . . (A15) 


The expression above gives the phase angle a of the radiated field (with 0 
and r fixed) as a function of the geometrical angle of observation 0 . In the case 
of a satellite spinning clockwise about the Z-axis 0 can be taken as* 


0 = 2 7T f t e Wt 

^ 8 S 


(A16) 


where f is the spin frequency of the satellite, and co g is the radian spin fre- 
quency. Then (15) becomes 


a = a Q + oj g t + a 1 sin(4o> | t + a 2 sin(8aj g t + + . 


(A17) 


^Taking the electromagnetic field of a rotating system to be identical to that of a fixed system 
implies the neglect of the field caused by mechanical motion of charges. This is permissible 
at low spin rates. 
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The one-way frequency shift resulting from this phase variation with time 
will be given by 



1 da 
2t 7 dt 


(A18) 


= f 8 + 4f s cos(4aj g t + yflj) 8f s a 2 cos(8a> s + /3 2 ) 


Equation (18) represents an addition to the Doppler shift already present 
because of the motion of the satellite in its orbit. The first term is a constant 
bias equal to the spin frequency, while the succeeding terms represent fluctua- 
tions at the 4th, 8th, etc. harmonic of the spin frequency. 

The absence of other harmonics is a direct result of the symmetry assumed. 
If, for example, the feeding currents l x and I 2 were not equal, then, referring 
to (A4) and .\(5), there would occur in place of (A9) and (A10) the results 


E„ (0 + 77 ) = - E a (0) = E e (0) el" (A19) 

E* (0 + 77) = - E* (0) = E 0 (0) el" (A20) 


and all even harmonics could appear in equation (A18). Complete lack of sym- 
metry, of course, would result in the appearance of odd harmonics. 

Simple examples of antenna patterns obtained with specific values of 
current distributions are given in the following sections. 


Pattern of a Small Turnstile Antenna 


A general expression for the field of a turnstile antenna is given in Appendix 
A, equations (A7) and (A8). If the antenna is small, i.e., if the length 4 of each 
monopole is small compared to the wavelength, then the currents on the mono- 
poles can be assumed to be approximately linear: 

f(x) = 1 -2L . (A21) 

4 
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Substituting this function into (A7) and (A8) in tho four integrals, approxi- 
mating each of the four cosine functions by one minus one half the square of its 
argument, and performing the indicated integration, there results, after some 
algebraic manipulation 


where 


and 


E, ? i - j B cos£ e J ^ [1 - 3b - b cos 40 + j b s in 40] 
E^, ? B [1 - b ■ b cos 40 - j b sin 40] 

B n Ij Tjki e-l kp / 4rrr 
b " k 2 i 2 sin 2 6/ 48. 


(A22) 

(A23) 

(A24) 

(A25) 


Since b << 1 for i << K , this can be put in the approximate polar form 


E 0 = B cos #(1 - 3b - b cos 4 0) e 


j (-’y +0+ b sin 40^ 


(A26) 


E* = B(1 - b + b cos 40) e '<*- b 8ln ''*) 


(A27) 


The phase angles of the components of the electric intensity agree with the 
general form (A15). 


Pattern of a Very Small Turnstile Antenna with Unbalance Feed 


The equations (A4) and (A5) give the field of a turnstile antenna fed by 
currents I 1 and I 2 . If the turnstile is very small, the currents on the mono- 
poles may be assumed, as in the section above to be given by equation (A21), 
and each integral in (A4) and (A5) is approximately equal to i/2 since each of 
the four cosine factors under the integral signs is approximately unity. With the 
further assumption that the feeding currents are not quite equal 



= q 2 1 


(A28) 
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but are in quadrature 


I^jqlj- (A29) 

The field of the antenna becomes approximately 

\ l fr**T llnM ) , Aom 

2B\e ' (A30) 

) - iJ L sin 2 $] 

e V 2 ' (A31) 

where B is given by (A24). 

As expected, in this case the phase pattern 


s B cos 6 ( 1 + q - 1 _ Izi cos 
0 1 2 2 


E. sBll + + lli cos 2 0 


J 4> 


2 2 


a(0) = 0 - 5 i. sin 2 0 (A32) 

2 

of exhibits a non-zero value of the second harmonic. 


CONCLUSIONS 

The general conclusion to be drawn from equation (A18) and from the ex- 
ample (A32) is that the one-way Doppler between the satellite and the ground 
station will be modified by the satellite spin with a constant bias equal to the 
spin rate* and by fluctuations at harmonics of the spin rate. Because of the 
electrical and physical symmetry of the antenna system, the fourth eighth, etc. 
harmonics should be substantial. The size oi the 2nd, 6th, etc. harmonics de- 
pends on how well the symmetry between the two antenna parts of the system is 
maintained in practice. Also, unless perfect symmetry exists, the fundamental 
and odd harmonics should be present to some small degree. 

As noted in the introduction, the satellites Explorer 33 and 35 actually em- 
ploy for tracking a two-way Doppler system in which transmission and reception 


*This possibility has been reported in the literature (Reference 11). 
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occur at different frequencies over the same antenna system. The system is 
balanced for transmission, and may be expected to exhibit some lack of balance 
for reception. The exact offect of such an arrangement on the amplitude of the 
bias and fluctuations depends on the frequency ratios and antenna connections 
used, and will be discussed in detail in a future report. The value of the bias 
will, however, be roughly the same as for the one-way case if the antenna con- 
nections are such that the spin frequency shift introduced in the uplink signal 
adds to the shift imparted to the downlink signal. For example,- if the spin rate 
f a is 1/2 cycle per second, then, assuming an approximate value of 2 meters 
for the VHF wavelength of the carriers, a range rate bias as large as (1/2) x2 = 
1 meter per second could oe obtained. The antenna connections used on Ex- 
plorers 33 and 35, however, are such that the frequency shifts caused by the- 
satellite spin in the uplink and downlink signal tend to cancel each other, and the 
net bias introduced is more than an order of magnitude smaller than the value 
of 1 meter per second given above. 

The frequency of fluctuations introduced by the spin, however, is the same 
as for one-way Doppler, no matter what connections are used. 
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APPENDIX B 


Method of Analysis 

For each data stretch, the following computations were performed. 

A least squares straight line T was fitted to the data, , 

T r a Q + a x i, (Bl) 

( 1 s 1 , 2, * • * n) 


and residuals, r i , from this straight line, 


r 


i 





were then analyzed for their spectral content.* 

Letting At be the time interval between adjacent values r A , mean lagged 
products c k , with lag interval- At =hAt, were calculated. 


n-hk 


q + hk 


■^5 V. 

q*Q 

^k = 0, 1, « • « , m, m < 


(B3) 


Raw spectral density estimates V k were obtained next. 


V k ,Ar 


* m- 1 

Cq + 2 y 

mtmmmml 

qd 


C q COS 


( qkrr\ 
m / 


+ C cos (kn) 


<B4) 


(k = 0, 1 m) 


*This particular method for obtaining the spectral density can be found in The Measurement of 
Power Spectra by R. B. Blackman and J, W. Tukey (reference 12). 



Finally, refined spectral density estimates U k wore calculated. 


U, = 1.08 V; + 0.92 Vj 


U. = 108 v m+ 0.92 

U k = 0.46 V, + 1.08 V k + 0.46 V k + l 


(B5) 


(k = 2, 3, * • * , m - 1) • 


Frequency Aliasing 


If a periodic function, say, cos (2?rf t),is sampled at equally spacer* intervals 
of time At, the following values of the function are obtained. 


g(i) s cos(2?r f i Ar) 


(i - • • • —2, — 1 » 0 » 1, 2i * * > ) > 


<B6) 


A sampling interval of A r seconds gives rise to a maximum resolvable 
(cutoff) frequency of f c <= (1/2 At) cycles/second. Letting f «= nf B be a harmonic 
of the satellite spin rate f # , g( i) can be written 


g(i) s cos(27m f 8 i At) 


5 cos £ < 2 77) (2) (k/2 Ar)i At ±27mf g iAT 3 


<B7) 


~ cos { ±2T7C(2kf c -nf s )iAr3) 

(k " 1, 2, * • * ) 

{i S * » * * Ij 0, 1 1 * * *) 

(rj ° 2, 3, * * ») 


n 



From B7 we see that for 1 2 k f c - nf s | < f c , the frequencies 1 2k f c - 
and nf g (k ~ 1 , 2, . . .) are indistinguishable and are called aliases of each 
other. 



APPENDIX C 


The Standard Deviation Of An Estimated Parameter 
I n A Least Squares Non-Linear Regression Model* 


Let z be an (nxl) vector of measurements of some variable such as two-way 
Doppler data, x a (kx 1) vector of parameters, and z (x) an(nx l) vector of com- 
puted values of the model which we assume adequately describes the observed 
data. We wish to determine that value x* of x which ilnimizes 


Q = (z - z(x) T <z - z(x)) 


(Cl) 


where the symbol T indicates the transpose, and k < n. ' 

For Q to be minimized, it is necessary that the gradient of Q with respect 
to x, evaluated at x = x*, vanish. Differentiating (Cl) 


F(x)=V x Q = -2A T (z-z(x)) (C21 

where 

/ 3zA 

A s j = ) = A (x) is an (n x k) mat rix of rank k < n 

\ dx J / 

By neglecting the dependence of A on x, we can use a modified Newton- 
Raphson iteration technique to find a solution x* such that F (x*) ~ 0. (See 
Figure C-l for a graphical representation of the Newton-Raphson technique 
in one dimension.) Differentiating (C2) gives 


F'(x) =7* (3 = 2A t A. (C3) 


The second derivative in (C3) always positive definite (since A was as- 
sumed to have rank k < n), Therefore, if the process cop marges to x% Q will 
be minimized, 


*TbU <Jov«>Ppm«nt con be found in Roforonco 13, 
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Figure Cl -Newton- Rapbson Iteration Technique In One Dimension 


If 


i 


Denote by x (m) the estimate of the solution at the m th Iteration. Then, the 
"improved" estimate x (m + is 

X<m+ 1 ) _ x (m) + A) ” 1 A T (z - Z (x (m) » (C4) 

where 

A ~ A(x< m >). 


Letting Ax = Ax (m + l) = (x (m+1 >- x (m >) and Az =Az (x< m >) = (z - z (x (m >)), 
we can rewrite (C4) as a "correction" to be applied at the (m + 1 ) th iteration. 


Ax = (A T A)” X A T Az 


(C5) 


where (A T A)~ X is the inverse of (A T A). 

Equation (C5) shows how small measurement errors (measured minus cal- 
culated) propagate through the assumed model in the neighborhood of the least 
squares solution x* = x (m + 1) . In other words, by neglecting the dependence of 
A on x, the error model indicated in (C5) becomes linear in the neighborhood of 

X* * x <m + l) # 


xt ..... 1 a. A rm La ** umuJ a uu fL /i ammam 4 xx f Kq wt aq oh norl 

JtNUW J.U1* {~\£j MU iX ruuuuiu VUUWl iupicugiiwiig VWO OJ-AVi U 1 mvuow*vM 


values z from the calculated least squares values z(x*) (the residuals) and AX 
a random vector representing the error in x due to AZ. We will assume that A 
has an expected value and covariance matrix given by 


Z 


8(AZ) = 0 

S( (AZ) (AZ) T ) =0-21 . 
where I is the ( n x n ) identity matrix. 


(C6) 


By (C6), we are assuming that the measurement errors are uncorrelated, 
have zero means and common variance p 2 , 


For the Goddard Range and Range Rate Tracking System, the sampling rates 
are less than the minimum bandwidth of the noise (10 Hz, reference 5). There- 
fore, if we have modeled the data adequately, having removed all deterministic 
trends, just "noise" with the properties in (C6) should remain. 
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Writing (C5) in random vector notation gives 


AX = (A T A)~ l A T AZ. 


<C7) 


From (C7) AX has expected value and covariance matrix, 

£(AX) = 0 

P = 8(AX A XT) =a 2 (A T A)- 1 . 


(C8) 


Since we do not know cr 2 , we must estimate it from the sample of residuals. 
Therefore, 


where 


and 


P= s 2 (A T A )’ 1 



(Z-Z(X*)) T (Z-Z(X*)) 


x* - x< m + 


<C9) 

(CIO) 


is an unbiased estimate of cr 2 in the neighborhood of X* where (by C5) the error 
model is linear (see reference 14 for the least squares estimator in the linear 
model). 

Therefore, the standard deviation <r xi . of X* (X* is one of the k components 
of X*) from (C9) is to a good approximation given by 



(Cll) 
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